Vortex stripe glass with self-generated randomness 
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We found a finite temperature glass transition in the absence of quenched disorder in frustrated 
Josephson junction arrays (JJA) on a square lattice under magnetic field with anisotropic Joseph- 
son couplings by numerical simulations. The externally induced vortexes develop zigzag stripes into 
the direction of weaker coupling at low temperatures. The whole amorphous vortex solid can flow 
smoothly into the direction of stronger coupling by non-trivial soft modes. As the result the macro- 
scopic phase coherence is destroyed being reminiscent of vortex flow in pure bulk superconductors 
or spin-chirality decoupling in frustrated magnets. On the contrary such soft-modes are absent into 
the direction of weaker coupling. Consequently the system behaves as a Ohmic (unjammed) liquid 
or superconducting (jammed) solid with respect to injection of external electric current, i. e. shear, 
along different directions. This jamming-glass transition can be regarded as a two dimensional 
realization of the Aubry's transition. 

PACS numbers: 61 .43.Fs,74.8f .Fa,74.25.Qt 



Glass transition is observed in diverse systems. An in- 
teresting basic question which continues to attract physi- 
cists is whether frustration of either energetic or kinetic 
origins, which are considered to be indispensable to avoid 
crystallisation, is sufficient to realize formation of amor- 
phous solids [l|, 2, 3]. The so called irrationally frustrated 
Josephson junction array (JJA) on a square lattice with 
irrational number density / of externally induced vor- 
texes is an ideal system to clarify this view [5| since it 
possesses strong energetic frustration in a simplest form. 
While JJAs with rational / develop periodic vortex lat- 
tices 0, 0: such simple orderings may be avoided for 
irrational /. Indeed equilibrium relaxations were simi- 
lar to the primary relaxations observed in typical fragile 
supercooled liquids [§]. 

Another remarkable feature of the JJA is that its elec- 
tric transport properties have formal analogy with the 
tribology and rheology : the current-voltage (IV) curves 
correspond to the flow curves, the shear-stress vs the 
shear-rate @. In a recent work we observed that 
variation of anisotropy A of the Josephson coupling into 
x and y-directions on the square lattice is analogous to 
the variation of the strength of interactions between in- 
commensurate surfaces in friction models II, IH that ex- 
hibit the so called Aubry's transition which is a jamming- 
transition at zero temperature. Such anisotropic JJAs 
of various A have already been fabricated by lithography 
techniques [l6| ■ In we have shown that the IV curves 
at T = with various A exhibit dynamical scaling fea- 
tures around the symmetric point A = 1 very similarly 



to the flow curves of granular systems [ll| around the so 
called jamming point (point- J) [l7j ]. 

In this letter, we present numerical evidences that the 
irrationally frustrated anisotropic JJA system exhibit a 
2nd order phase transition into a glassy phase at a finite 



temperature T C (A). Since there is no quenched disorder 
in the present system, the glass transition here is fun- 
damentally different from those in the spin-glasses and 
conventional vortex glasses [6J but closer to structural 
glass transitions of liquids. T C (A) rapidly decreases with 
the decrease of the anisotropy |A — 1| and quite possi- 
bly terminates at the 'point-J' (T = 0, A = 1), being 



consistent with recent studies [14|, [15( which suggest the 
isotropic system with A = 1 remains in the vortex liquid 
state at any finite temperatures T. 

Model The JJA [j, [H is a network of super- 
conducting islands. We consider a square lattice of size 
N = L 2 and label the islands as i = 1, . . . , N. The phase 
of superconducting order parameter 0i of the i-th island 
is coupled to its nearest neighbors by the Josephson junc- 
tions. External transverse magnetic field B threads the 
cells in the forms of flux lines each of which carrying a 
flux quantum (f>o — ch/2e. On average each unit cell with 
area a 2 carries / = Ba 2 /<fio flux lines. A flux line thread- 
ing a unit cell induces a vortex of the phases around the 
cell much as a dislocation in a crystal. The number den- 
sity / is chosen to be an irrational number. Thus the 
system is extremely frustrated such that even formation 
of periodic vortex (dislocation) lattice can be avoided. 
The properties of the JJA under the transverse magnetic 
field are known to be described by an effective Hamilto- 
nian [HH, 
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cos(V>tf)-A cos(V>ij) (1) 

<ij> || x-axis 



<ij>||j/-axis 



with the gauge-invariant phase difference, Vij = &i — &j — 
Aij. We measure temperature T in a unit with = I. 
Here A is the ratio of the strength of the Josephson cou- 
plings along x and y-axis. In the following we only ex- 
amine A > 1, which is obviously sufficient by symmetry. 
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FIG. 1: Vortex patterns and the structure factor. Here 
A = 1.5 and T = 0.2 (T c ~ 0.23). a) display the vortex 
pattern in a snapshot in equilibrium and b) displays that at 
a nearby energy minimum reached via an energy descent al- 
gorithm, c) displays the cross-sections of the structure factor 
S(q x ,q y ) (L = 55) with thermal average (thick lines) and at 
the energy minimum shown in a) (thin lines) . d) displays the 
amplitude of the peak of the structure factor S* = S(n, 2nf) 
with thermal average (S»), average over the energy minima 
[S^lmm. and variance of the minimum-to-minimum fluctuation 
V[S?]min. - [S*]min.- Here [. . .] min is took over 100 energy- 
minima obtained by independent initial conditions. 

The vector potential A»j(= —Aji) is defined such that 
directed sum of them around each cell is 2irf. 

Charge Vi of the vortex at the cell associated with the i- 
th lattice site is defined by taking directed sum of (ipij — 
[tpij] n ) /2tt on the junctions around the cell. Here [x] n 
denotes the nearest integer of the real variable x. It takes 

values . . . , — 1 + /, /, 1 + /, The total charge is zero 

due to the charge neutrality. 

Methods In numerical simulations, we used a series of 
rational numbers p/q = 5/13, 8/21, 13/34, 21/55, 34/89, 
55/144, 89/233 for the filling /, which approximate an 
irrational number / = (3 - a/5)/2=0.38196601... We 
worked on systems of size L x L with L = q so that in 
the thermodynamic limit L — ► oo the value of / converges 
to the target irrational number. 

To investigate macroscopic static properties, we used 
Monte Carlo (MC) simulations. We performed simula- 
tions on L = 13 — 89 using 20 — 120 temperatures in the 
temperature range T = 0.2 — 0.4 by the exchange MC 
method [HI]. More precisely each MC step consists of 
one sweep over the system by the Metropolis updates of 
each phase 0,, one step of the overrelaxation [24j and one 
step of the exchange MC. We used 10 5 - 10 6 MC steps 



FIG. 2: Correlation lengths and two-time auto-correlation 
functions, a) Correlation length of the vortex £v and c) phase 
£p. Here A = 1.75 for which T c ~ 0.25. b) Auto-correlation 
functions of vortex Cy(t + i w , t w ) and d) phase Cp(i + i w ,£ w ) 
during aging at T = 0.2 after temperature quench with wait- 
ing times t w = 10,10 2 ,10 3 ,10 4 ,10 5 and 10 6 from left to right 
curves. 



for the equilibration and observations. 

We also measured two-time auto-correlation functions 
during relaxation from random initial configurations by 
a simple Metropolis updates over 10 6 MC steps on larger 
systems with L = 144, 233 which do not show appreciable 
finite size effects within the simulation time scale. 

Vortex stripe glass As shown in Fig. 1 a), the vor- 
texes develop zigzag folded layers along the direction of 
weaker coupling which are stacked along the direction of 
stronger coupling at low temperatures. The formation 
of the smectic like stripes is reasonable because the re- 
pulsive interactions between vortexes are anisotropic. A 
remarkable feature is that the zigzag folds at different 
layers are strongly correlated. Indeed, as shown in Fig. 1 

b) , the alignment of the zigzag pattern is nearly perfect 
in the nearby energy minimum which we obtained by 
a simple energy descent algorithm. As shown in Fig. 1 

c) d), the structure factor S{q Xl q y ) exhibits prominent 
peaks at (q x ,q y ) — (n,2irf) and (tt, 2tt(1 — /)) whose 
height scales with the system size as N = I? as usual 
Bragg peaks. However the profile of the peak is pecu- 
liar: it decays sharply along q y but decays slowly by a 
power law \ir— q x \~ 2 along q x . The exponent is essentially 
independent of the temperatures. These features reflect 
the random zigzag wandering of layers which are strongly 
correlated between different layers. Somewhat different 
smectic layering of vortexes exist in bulk high-Tc super- 
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conductors due to the presence of the oxide layers [2l| . 

Fig. 2 a) displays the scaled correlation length £.y /L 
of the vortexes measured by the same method as in [15| . 
(See [l5[ for the details of the method, which is a stan- 
dard technique for spin-glasses.) Here the correlation 
length along the direction of stronger coupling is shown. 
The scaled correlation length £v /L of different sizes ex- 
hibit crossings suggesting a 2nd order phase transition at 
TJX), at marked variance with the isotropic A = 1 case 
[H in which £ v ~ T- 2 - 2 . 

Absence of phase long-ranged order For the case 
of isotropic coupling (A = 1), it has been noticed that 
correlation length of the phase grows less rapidly £p ~ 
T~ 12 than that of the vortex £ v ~ T~ 2 ' 2 . This is very 
similar to the spin-chirality decoupling observed in a 2 
dimensional XY spin-glass [23[ in the sense that chirality 
is equivalent to the externally induced vortex. 

We found the decoupling is much more enhanced in 
the presence of anisotropy. Fig. 2 b) displays the scaled 
correlation length £p/L along the direction of stronger 
coupling. The data do not exhibit crossings at larger L. 
Furthermore the two-time phase auto-correlation func- 
tion C P (t+t w ,t w ) = (l/N) E,=i( cos (0<(<+*w)-0<(*w))> 
do not exhibit any plateau as shown in Fig. 2 d). These 
features suggest absence of phase long-ranged order. 

Ergodicity breaking Starting from different initial 
conditions we could obtain numerous energy minima sim- 
ilar to the one shown in Fig. 1 b) but with different ran- 
dom configurations of the zigzag fold. The randomness is 
reflected for example in the minimum-to-minimum fluc- 
tuations of the structure factor shown in Fig. 1 c) d). 
Note that both the average and the variance of the fluc- 
tuation grows linearly with the system size N. Quite re- 
markably we found the energy differences among them 
are only O(l) so that all of them are relevant in the static 
equilibrium ensemble. However, the nearly perfect ver- 
tical alignment of the zigzag fold suggests that it would 
take extremely long time to change the pattern of the 
zigzag fold from one to another dynamically. 

To clarify this point, we performed aging simulations 
and measured the two-time vortex auto-correlation func- 
tion CyQt + tw, *w) = (l/#)££i<«i(* + **)«<(*))• As 
shown in Fig. 2 b) it exhibits a plateau at larger waiting 
times t w , which suggests the Edwards- Anderson order 
parameter of the vortex q\? = lim^oo lim^^oo Cy(t + 
t w ,t w ) is finite, i. e. the ergodicity is broken. It implies 
that once the system reaches an energy minimum, as the 
one shown in Fig. 1, it cannot escape out of it. Thus in 
sharp contrast to the usual smectic liquid crystals [20| . 
the transverse wandering of the layers are frozen in time. 
Note that usual smectic liquid crystals have only quasi- 
long raged orders along all directions so that the layers 
are not so regularly stacked as in the present system. 

Soft modes To obtain some insights on the low lying 
metastable states, we analyzed the spatial configuration 
of the Josephson currents sin(^jj). Since the vortexes 
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FIG. 3: Spatial configuration of currents through the Joseph- 
son junctions in an energy minimum (shown in Fig. 1 b)). Co- 
ordinate along either x or y axis is "folded" to elucidate mod- 
ulation with respect to ideal incommensurate ordered states 
(see text). Panels a) b) show Josephson currents into x and y 
directions plotted against the folded y-axis. Similarly panels 
c) d) show the currents against the folded :r-axis. 



have long ranged order with the incommensurate wave 
number q y = 2nf into the direction of stronger coupling, 
which is the y-direction for A > 1, we may take ideal 
incommensurate ordered states as references. Then it 
is convenient to introduce a "folded coordinate" fy — 
Int(fy) with Int(x) being a function which truncates a 
real number x to an integer number [26j |. 

Fig. 3 a) b) display the configuration of the current J x 
and J y along the junctions parallel to the x and y axises. 
Apparently the currents behave as a continuous function 
of the folded y-axis. This immediately means existence of 
soft-modes: by choosing an arbitrary real number a and 
changing the entire configuration such that the value of 
the gauge invariant phase difference ip at (x, fy— Int(/y)) 
is replaced by that at (x, f(y+a)—hxt(f(y+a)) we obtain 
a different configuration. By a little thought one can find 
that it is another energy minimum with exactly the same 
energy. We note that J x takes all possible values between 
— 1 and 1 while J y takes values only within a narrow range 
around 0. It means that the shift by a certain a in the 
folded y-coordinate amount to put shear across the x- 
direction by which the whole vortex solid slide along the 
y-direction without changing its pattern and the energy: 
in spite of the fact that the vortexes live on the lattice, 
they flow like in the pure bulk superconductors fl8j |. 

In contrast, the configuration of the currents with the 
folded x-axis do not reveal any continuous functions but 
discontinuous ones as shown in Fig. 3 c) d). This means 
the vortex solid cannot slide freely along the ir-direction. 

The presence/absence of the soft modes into differ- 
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FIG. 4: Phase diagram and helicity modulus, a) Note that 
T c (l/A) = T C (X)/X holds due to the obvious symmetry be- 
tween x and y axis, b) The profile of helicity modulus Y x and 
Y y for shear along x and y directions at A = 1.75. 



surement of the relaxation times Ty(T, A) of the vor- 
tex auto-correlation function Cy(t) which can be fit- 
ted as r v (T,A) ~ (T - T c {\))- Zl/ with zv ~ 6.5 and 
T C (A) ~ |A- l| p with p ~ 0.4. More details of the results 
will be reported elsewhere. 
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ent directions naturally explain the unusual features ob- 
served in [l3| that shear (helicity) modulus is finite/zero 
with respect to shear (injection of external electric cur- 
rent) along the direction of stronger/ weaker coupling 
(See Fig. 4 b)). The isotropic point A = 1 appears 
as a critical point ('J-point') for the shear response at 
T = 0. Correspondingly dynamic shear response also 
exhibit jammed/unjammed behaviours as shown in [loj ]. 

Furthermore the phase coherence of the system along 
the direction of weaker coupling must be destroyed by 
the soft mode by which the vortex solid slide into the 
direction of the stronger coupling. Thus the phase co- 
herence can only grow in one dimensional way along the 
direction of stronger coupling resulting in the absence of 
long-ranged phase order, as observed here. 

Aubry's transition as a jamming-glass transi- 
tion In an one dimensional limit the 2 dimensional JJA 
reduces to a JJA on a two-leg ladder lattice [22j], which 
can be viewed as an Eulerian representation of models 
for friction between two atomic layers 12, IH . The ra- 
tio of the atomic spacings of the layers is assumed to be 
an irrational number /, which is identified to the vortex 
density / in the JJA after some gauge transformations. 

The Josephson currents sin(^>y) shown in Fig, 3 cor- 
respond to inter/intra layer forces between atoms in the 
friction models. In the Frenkel-Kontorova model it is 
known rigorously that they are analytic functions of 
the folded coordinate in the unjammed phase but non- 
analytic functions with dense set of discontinuous points 
in the jammed phase where the static frictional force or 
yield stress (corresponding to the critical current in JJA) 
is finite (T^. They are called as the Hull functions. 

Our results discussed above strongly suggest the 
jamming-glass transition we found in the present model 
is a two-dimensional analogue of the Aubry's transi- 
tion. In contrast to the one-dimensional models, how- 
ever, the glassy phase is observed not only at zero tem- 
perature but also at finite temperatures. In Fig. 4 a) 
we display the phase diagram of the present system. 
The critical temperature T C (X) is estimated by mea- 
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